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Abstract. Several new aspects of the subtle interplay between electronic correla- 
tions and disorder are reviewed. First, the dynamical mean-field theory (DMFT) 
together with the geometrically averaged ("typical") local density of states is 
employed to compute the ground state phase diagram of the Anderson-Hubbard 
model at half-filling. This non-perturbative approach is sensitive to Anderson lo- 
calization on the one-particle level and hence can detect correlated metallic, Mott 
insulating and Anderson insulating phases and can also describe the competition 
between Anderson localization and antiferromagnetism. Second, we investigate 
the effect of binary alloy disorder on ferromagnetism in materials with /-electrons 
described by the periodic Anderson model. A drastic enhancement of the Curie 
temperature T c caused by an increase of the local /-moments in the presence of 
disordered conduction electrons is discovered and explained. 



1 Introduction 

Electronic correlations are known to lead to a plethora of fascinating phenomena [1I2I3I4I5I6I7I8] . 
The same holds true for disorder in quantum systems [911011 1112113114115116117] . Both limits, 
electronic correlations without disorder and disorder without electronic correlations, are notori- 
ously difficult to treat since interactions lead to a highly complicated many-body problem, while 
disorder requires the application of statistical methods for taking averages. The simultaneous 
presence of interactions and disorder, often encountered in real materials, therefore implies an 
even more complex many-body problem [18I19I20I21I22I23I24I25I26I27I28I29"I3"0] which is still 
far from understood. Namely, repulsive interactions and disorder are both driving forces be- 
hind metal-insulator transitions (MITs) connected with the localization and delocalization of 
particles. While the Coulomb interaction may trigger a Mott- Hubbard MIT |4I6I31| . the scat- 
tering of non-interacting particles from randomly distributed impurities can lead to Anderson 
localization [9132] . 

The exploration of phenomena which take place at intermediate or strong couplings, e.g., fer- 
romagnetism and the Mott-Hubbard MIT with and without disorder are of particular interest. 
Here the recently developed dynamical mean-field theory (DMFT) |33I34I35I36I37I38I39I40I4T] 
has proved to provide an excellent, comprehensive mean-field approximation, which may be 
employed for arbitrary values of the input parameters. For this reason the DMFT has been 
successfully applied in the investigation of electronic correlation effects in theoretical models 
and even real materials [41)42143144] . 



a e-mail: dieter.vollhardt@physik.uni-augsburg.de 



2 



Will be inserted by the editor 



In this paper we review our recent results on correlated electron systems in the presence 
of disorder obtained within the DMFT. In particular, we explore the low-temperature phase 
diagram of the Anderson-Hubbard model at half-filling, where disordered metallic, Anderson 
localized, and Mott insulating paramagnetic phases compete with antiferromagnetic long-range 
order. For comparison we then investigate the phase diagram of the simpler Falicov-Kimball 
model. Finally we discuss the influence of the band-splitting by binary alloy type disorder on 
properties of the periodic Anderson model. In particular, we show that alloy disorder not only 
leads to a non-monotonic changes of the Curie temperature T c as a function of some control 
parameter, and even to an enhancement of T c compared to the non-disordered case, but also 
to the formation of Mott or Kondo insulators at non-integer electron densities. 



2 Correlated electron systems in the presence of disorder 

The models investigated here are the Anderson-Hubbard model 

H =Y1 *« c L c j CT + E eiHia + U Y n ^ n ^ (!) 

ija ia i 

where is the hopping matrix element, U is the local Coulomb interaction, c\ a (ci a ) is the 
fermionic creation (annihilation) operator for an electron with spin a in Wannier state i, and 
nj CT is the particle number operator; and the Anderson-Falicov-Kimball model 

h = tv c k + E e ^ + uJ2 fUA^ (2) 

ij i i 

where cf (//) and Cj (/J are fermionic creation and annihilation operators for mobile (immobile) 
particles at a lattice site i. Furthermore, is the hopping amplitude for mobile particles be- 
tween sites i and j, and U is the local interaction energy between mobile and immobile particles 
occupying the same site. The ionic energy ej in both models is a random, independent variable 
which describes the local, quenched disorder affecting the motion of the mobile particles. 

The disorder part is modeled by a corresponding probability distribution function (PDF) 
P(e^). A model of disorder that we use for studying Anderson localization is one with the 
continuous PDF 

with (9 as the step function. The parameter A is then a measure of the disorder strength. 



2.1 Dynamical mean-field theory with disorder 

Dynamical mean-field theory is based on the observation that in the high-dimensional limit 
d — > oo (or equivalently Z — > oo, where Z is the lattice coordination number) the self-energy 
Uij(uj) as defined by the Dyson equation 

GjjffK) -1 = G^vK)- 1 - (4) 

where i,j denote lattice sites and w n = (2n + 1)tt//3 are fermionic Matsubara frequencies, 
becomes diagonal in real-space |33l34j 

£ijcr(Wn) = Eicr(Un) (5) 

provided the hopping amplitudes are scaled properly to ensure a balance of kinetic and interac- 
tion energy in this limit. For nearest-neighbour hopping on a hypercubic lattice this implies the 
scaling t = t*/V2d = t* j\[Z. In a homogeneous system the self-energy is also site independent, 
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i.e., Sij a (iu} n ) = E a {iui n ) Sij, and is only a function of the energy. The DMFT approxima- 
tion when applied to finite-dimensional systems neglects off-diagonal parts of the self-energy. 
In other words, the DMFT takes into account all temporal fluctuations but neglects spatial 
fluctuations between different lattice sites [41143] . 

Here we present recent developments regarding the application of DMFT to correlated 
fermion systems with disorder. Within DMFT each correlated lattice site is mapped onto a 
single impurity, which is coupled to a dynamical mean-field bath describing the influence of all 
remaining lattice sites. This coupling is represented by the hybridization function i]i a (oj), which 
needs to be determined self-consistently. The mapping is performed for all Nl lattice sites. 

In the presence of disorder as given by a particular realization of onsite energies {ei, €2, £n l } 
the total partition function can be written as a product of the Nl partition functions of the 
individual impurities 



z= n z ' = n cx p(E in [* 

l i \(7LO n 



Via(u n ) - E ia (u> n )} . (6) 



where the hybridization function rji a (ijj n ) formally represents a site- and time-dependent one- 
particle potential. The unitary time evolution due to this potential can therefore be described 
by a local, time-dependent evolution operator [45146] 

where the interaction representation has been used, T T is the time ordering operator, and Ci CT (r) 
evolves according to the atomic part H ' oc of the Hamiltonians ([1]) or © in imaginary Matsubara 
time r G (0,/3). We write the partition function (0 as a trace 

N L 



Z = Z[r Ha ] = [] Tr [e^^-^^]] , (8) 



where Nj oc is the local particle number operator. 

For a given dynamical mean- field rji a (Lo n ) the local one-particle Green function Gu^^uin) 
is then determined from Eq. (|5|) by taking a functional logarithmic derivative of the partition 
function ([5]) with respect to rj^ (uj n ) 

We thus obtain the local Dyson equations 

Si„((J n ) =lUJ n + fl- €i - Tfoiun) - 7T—( — X > (10) 

for each of the Nl lattice sites. Eqs. (jH [3 13 HI and lT0|) constitute a closed set of self-consistency 
relations. The solution of these yields an approximate solution of the Hamiltonian ([1]) or 
for a given disorder realization. 



2.2 Disorder averages 

Solving Eqs. (HI [5j HI [9j and ITUl) rigorously is in general only possible for small Nl, since it 
requires an exact evaluation of the time evolution operator ([7]) at each impurity site. However, 
in disordered systems one is mostly interested in the thermodynamic limit Nl — > 00 where 
localized and extendend states can be distinguished reliably. Here one faces a typical trade- 
off situation in computational physics. A solution to this problem is obtained by a statistical 
approach in combination with DMFT, as applied by us [47148149] and outlined in the following. 
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Given a particular disorder realization {ei, 62, £n l } one could in principle obtain from 
the solution of the DMFT equations a set of realization-dependent local densities of states 
(LDOS) 

A ia {u) = ImG !W (u) ffi -> ui + i0 + ) . (11) 

Usually, however, one is interested in physical information about a system that does not depend 
on a particular disorder realization. This makes a statistical interpretation of the solutions of 
Eqs. (gl E El and EH! necessary. 

A common approach for systems in the thermodynamic limit Nl ~ > 00 is to take the 
arithmetic average of the LDOS Ai a {uS) over many realizations of the disorder, i.e., 




By performing the arithmetic average one restores the translational invariance in the description 
of the disordered system, i.e. the average i4 er (o;) ar ith = (Ai a (u)) is the same for all lattice sites. 
However, this only leads to meaningful results if both the system and the physical observable 
under study are self-averaging. 

An example of a non-self-averaging system is a disordered system at the Anderson localiza- 
tion transition |9 j . This implies that during the time evolution, a particle cannot explore the full 
phase space, i.e., cannot probe all possible random distributions. In such a case the arithmetic 
average (fT2"]) is inadequate. Here one is faced with the question concerning the proper statistical 
description of such a system. 

As pointed out by Anderson [5] , the solution to this problem is to investigate the full PDF 
for a given physical observable P[Ai a (u>)] and to identify its most probable value - the "typ- 
ical" value A CT (a;)typ - i.e. the value where the PDF P[Aj CT (a;)] has a global maximum. This 
typical value will be the same on each lattice site and represents typical properties of the sys- 
tem. By considering A CT (w)ty P one thus restores translational invariance in the description of 
a disordered system. Using photoemission spectroscopy one could, in principle, experimentally 
probe the LDOS at different lattice sites and measure its most probable value. We note that 
if sample-to-sample fluctuations are small, the typical value ^4 CT (o;)typ coincides with the arith- 
metic average Ar(w) a rith- On ^ ne other hand, in a non-self-averaging system the PDF can be 
strongly asymmetric, with a long tail, in which case the typical value A a (ui) tTP would strongly 
differ from A(j(cj) ar ith- In this case the arithmetic mean is strongly biased by rare fluctuations 
and hence does not represent the typical property of such a system. 

Calculating the full probability distribution of the LDOS typically requires the inclusion of a 
large number of impurity sites, which for interacting systems is hard to achieve computationally, 
although there have been recent successful attempts in this direction [50] ■ A more efficient 
approach is based on identifying a generalized average which yields the best approximation to 
the typical value. Among several different means the geometric mean 

AxMgcom = exp [(In A ia {u))} , (13) 

turns out to be very convenient to describe Anderson localization, where it represents a good 
approximation to the typical value: 

ArMtyp « A a (uj) gcom . (14) 

Here (-F(ej)) — J Yl i deiP(ei)F(ei) denotes the arithmetic mean of the function F(ei). It is 
easy to see that in the special case when P[Ai a (uj)\ is a log-normal distribution the relation 
A IJ {uj)t Y p — A a (ijj) geom holds exactly. For the noninteracting Anderson disorder model it was 
shown that A CT (u;)geom vanishes at a critical strength of the disorder, thereby providing an 
explicit criterion for Anderson localization [9I51I52I53I54J . 

By using the geometrically averaged LDOS one arrives at a translationally invariant de- 
scription of the disordered system. This allows one to solve the DMFT equations in the ther- 
modynamic limit as will be demonstrated in the next section. The problem of finite-size effects 
is then absent, and the main limitation is due to the finite accuracy of the impurity solver. 
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2.3 Self-consistency conditions of the DMFT with disorder 

Once the geometrically averaged local spectrum has been obtained, the corresponding local 
Green function is given by 

G CT (a, n ) geom = / a /^W , (15) 

J lLO n - LO 

which allows one to recast the DMFT self-consistency condition ([TU]) into the following trans- 
lationally invariant form 

S a {ui n ) = iui n + fi - ri a (uj n ) - 1 . (16) 

tr y^n J geom 

where a vanishing average on-site energy (ej) = has been assumed, which holds in particular 
for the box-shape PDF. On the other hand, Fourier transformation of the lattice Dyson equation 
© yields 

G CT (w n ) g eom = / dz- r, (17) 

which closes the set of DMFT self-consistency relations. We note that the self-consistency 
equation (|17[) gives the same geometrically averaged Green function as that obtained from 
([15)) . Here Nq(z) is the density of states for a non-interacting lattice system without disorder. 

The above discussion shows that by taking the geometric average of the LDOS one i) restores 
translational invariance and ii) can directly address the thermodynamic limit. The DMFT 
calculations are thus no longer affected by finite-size effects. 

For the description of antiferromagnetic long-range order the self-consistency relations need 
to be modified. In this case one introduces two sublattices s =A or B, and calculates the 
two corresponding local Green functions Gii as (uj n ), which are no longer identical. Geometric 
averaging of the LDOS given by Eq. (fTTl) yields ^4 CTS (a;)gcom = Gxp[(\n A icrs (uj))}. The local 
Green function is then again obtained from the Hilbert transform (I15p . and the self-energy 
£ as (uj) from the local Dyson equation Eq. (TIB)) . Finally, the self-consistent DMFT equations 
are closed by the Hilbert transform of the Green function on a bipartite lattice: 

G as (u) n ) S eom = I dz -= N o{ z ) _ _ 

where s denotes the sublattice opposite to s [55 41] . 

As a concluding remark regarding the formalism we note that replacing the geometric mean 
by arithmetic averaging would lead to a description of disorder effects on the CPA level only, 
which would not allow to detect Anderson localization. We also point out that in the presence 
of disorder the typical LDOS A a (uj) gcom is not normalized to unity. This means that ^(ccOgeom 
only describes extended states, i.e. the continuum part of the spectrum. The contribution of 
localized states to the LDOS is not captured by DMFT with geometric averaging. Therefore, 
this approach cannot describe spectral properties of the Anderson-insulator phase. 



3 Phase transitions in the disordered Hubbard model 

3.1 Characterization of phases 

Our goal is to determine the ground-state phase diagram of the Anderson-Hubbard Hamilto- 
nian (fT]). In order to characterize the different relevant phases, we compute several physical 
observables: 

1. the LDOS ^4(Ts(w)geom for a given sublattice s and spin direction a; 

2. the total DOS for a given sublattice s at the Fermi level (uj = 0) with N s (0) geom = 
J2<, A <rs{u = 0) g oom; 
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3. the staggered magnetization m^™ 1 = \n^ m — n^™ m |, where n|° om = J ^ dwA as (uj) gcoln is 
the local particle density on sublattice s. 

The possible phases of the Anderson-Hubbard model can then be classified as follows: The 
systems is a 

— paramagnetic metal if Nf eom (0) ^ and rn^p" 1 = 0; 

- AF metal if Af com (0) ^ and mf° m ^ 0; 

- AF insulator if jVf eom (0) = and mfp™ ^ but Nf eom (u) ± for some w ^ (in fact, the 
last condition is already implied by m^ C p ln ^ 0); 

— paramagnetic Anderson-Mott insulator if Nf eom (w) — for all u>. 

In the following we consider the Anderson-Hubbard model at half-filling and determine 
its paramagnetic as well as magnetic phases at T = using the formalism described above 
|47l49j . For computational convenience, we choose a model DOS, N (e) = 2y/D 2 — e 2 /nD 2 , 
with bandwidth W — 2D, and set W = 1 in the following. For this DOS and for a bipartite 
lattice the local Green function and the hybridization function are connected by the simple 
algebraic relation 77 <TS (w) ge om = -D 2 G cr s(w) g0 om/4. gl] 

The resulting DMFT equations are solved at zero temperature by the numerical renormal- 
ization group technique [56 57 58 , which allows us to solve the effective Anderson impurity 
problem and thus to determine the LDOS and its averages. 

3.2 Mott-Hubbard transition vs. Anderson localization 

In Fig. Q] we show the paramagnetic ground-state phase diagram of the Anderson-Hubbard 
model at half-filling obtained by the DMFT formalism detailed above. Two different quantum 
phase transitions can be observed: a Mott-Hubbard MIT for weak disorder A, and an Anderson 
MIT for weak interaction U. In the following we will discuss these transitions as well as the 
properties of the two different insulating phases and of the correlated, disordered metallic phase. 
Since these phases are all paramagnetic the spin index will be omitted in this section. 

(i) Disordered, metallic phase: The correlated, disordered metal is characterized by a finite 
value of the spectral density at the Fermi level, A(u> = 0) gO om 7^ 0. 

(ii) Mott-Hubbard MIT: For weak to intermediate disorder strength a sharp transition be- 
tween a correlated metal and a gapped Mott insulator occurs at a critical value of U. Two 
transition lines are found depending on whether the MIT is approached from the metallic 
[A^ H (U), full dots in Fig. [T] or from the insulating side [A^( H (U), open dots in Fig. [J. The 
curves A^{ H {U) and A^ H (U) in Fig. Q] are seen to have positive slope. This is due to the 
disorder-induced increase of spectral weight at the Fermi level which in turn requires a stronger 
interaction to open the correlation gap. As a result, close to the hysteretic region an increase 
of disorder will drive the system from the Mott insulator back into the metallic phase or, 
equivalently, protect the metal from becoming a Mott insulator. 

The transition lines A^[ H (U) and A^ H (U) terminate at a single critical point around 
A fa 1.8, cf. Fig. [1] For stronger disorder (A > 1.8) we observe a smooth crossover rather than 
a sharp metal-insulator transition. 

(Hi) Anderson MIT: Next to the metallic phase and the crossover regime an Anderson in- 
sulator is found where the LDOS of the extended states, as obtained by geometric averaging, 
vanishes completely (see Fig. [TJ . The associated critical disorder strength Af(U) at the Ander- 
son MIT has a non-monotonic behaviour: while it first increases as a function of interaction in 
the metallic regime, it decreases again in the crossover regime. In the former case, an increase 
of interaction strength leads to re-entrant transition from the Anderson insulator into the cor- 
related metal, where electronic correlations prevent quasiparticles from localization by elastic 
impurity scattering. 

(iv) Anderson-Mott insulator: We note that the Mott insulator (with a finite correlation 
gap in the single-particle spectrum) is rigorously defined only in the pure system (A = 0), while 
the gapless Anderson insulator phase arises only for non-interacting systems (U = 0) with suf- 
ficiently strong disorder A > A^(0). In the simultaneous presence of finite interactions and 
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Fig. 1. Paramagnetic ground-state phase diagram of the Anderson-Hubbard model at half-filling 
obtained by DMFT with the typical local density of states 47 . 

disorder it is no longer possible to strictly distinguish these two phases. However, as long as the 
LDOS shows the characteristic Hubbard subbands the system may be termed a disordered Mott 
insulator. With increasing disorder A the spectral weight of the Hubbard subbands vanishes 
and the system can be considered a correlated Anderson insulator. The (crossover) boundary 
between these two types of insulators is marked by a dashed line in Fig. [TJ Our DMFT results 
obtained here indeed show that the paramagnetic Mott and Anderson insulators are continu- 
ously connected. Hence, by changing U and A it is possible to move from one insulating state 
to another one without entering a metallic phase. This single connected, insulating phase is 
therefore termed the Anderson-Mott insulator. 

3.3 Competition between Anderson localization and antiferromagnetism 

In the previous section we have neglected antiferromagnetic ordering, which generically occurs at 
low temperatures on non-frustrated lattices and in the absence of disorder. Several fundamental 
questions arise: (i) How do local interactions influence a non-interacting, Anderson localized 
system at half-filling? (ii) How does an antiferromagnetic (AF) insulator at half-filling respond 
to disorder which in the absence of interactions would lead to an Anderson localized state? (iii) 
Can Slater and Heisenberg antiferromagnets be distinguished by their response to disorder? 
Here we provide answers to these questions by calculating the zero temperature, magnetic 
phase diagram of the disordered Hubbard model at half-filling using DMFT together with a 
geometric average over the disorder and allowing for a spin-dependence of the DOS [49] . 

Our results are shown in Fig. [5] The response of the system to disorder is found to be 
qualitatively different depending on whether the interaction U is weak or strong. At strong 
interactions, U/W > 1, only two phases exist, an AF insulating phase at weak disorder, A/W < 
2.5, and a paramagnetic Anderson-Mott insulator at strong disorder, A/W > 2.5. The local 
DOS and the staggered magnetization both decrease gradually as the disorder A increases 
and vanish at the phase boundary, indicating that the associated quantum phase transition 
is continuous. On the other hand, a richer structure of the phase diagram is found for weak 
interactions, U/W < 1, (Fig. 0). In particular, for weak disorder a paramagnetic metallic phase 



8 



Will be inserted by the editor 



3 



PM Anderson-Mott insulator 




2.5- 



AF insulator 



Heisenberg gap 



0, 







0.5 



1 1.5 

Interaction U 



2 



2.5 



Fig. 2. Magnetic ground-state phase diagram of the Anderson-Hubbard model at half- filling as obtained 
by DMFT with a spin resolved, geometrically averaged local DOS (see text) [49] . 

is stable. It is separated from the AF insulator at large U by a narrow region of an AF metallic 
phase. The AF metal is characterized by long-range order in the absence of a gap, which is due 
to a redistribution of spectral weight induced by disorder [49] . 

In the pure system without disorder, the AF insulating phase can be characterized by two 
limiting regimes: a Slater antiferromagnet with a small gap at U /W < 1 and a Heisenberg regime 
with large gap at U/W > 1. These limits can be addressed by perturbation expansions in U 
and 1/U around the symmetry broken state of the Hubbard and the corresponding Heisenberg 
model, respectively. Our results for toaf confirm that even in the presence of disorder there 
is no sharp transition between these limits, in agreement with earlier studies [59) . This may 
be attributed to the fact that both the Slater and the Heisenberg regime are characterized by 
the same order parameter. However, from our results (Fig. [2]) it is evident that the two limits 
can be distinguished by their overall response to disorder. Namely, the reentrance of metallic 
antiferromagnetism at A/W > 1 occurs only within the Slater AF insulating phase. 

Let us finally remark that a disordered AF metallic phase is also obtained by DMFT com- 
bined with arithmetic averaging [55160] . This approach however predicts that both the param- 
agnetic and the AF metal remain stable for arbitrarily strong disorder, which is clearly incorrect 
and closely related to the failure of CPA to capture disorder-induced localization. Only a com- 
putational approach which is sensitive to Anderson localization, such as the DMFT with geo- 
metrically averaged local DOS employed here, is able to detect already on the one-particle level 
the suppression of the metallic phase for A/W > 1.36 and the appearance of the paramagnetic 
Anderson-Mott insulator at large disorder A. 

4 Metal-insulator transitions in the disordered Falicov- Kimball model 

We now apply the DMFT with geometric average over the disorder to the disordered Falicov- 
Kimball model at half- filling [?5]. The non-interacting DOS for mobile particles and other 
physical parameters are the same as in the previous section. However, the DMFT equations 
are now solved differently: since in this case the functional integrals are Gaussian they are 
calculated analytically, and only the self-consistency loops are computed numerically. 
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Fig. 3. Ground-state phase diagram of the Falicov-Kimball model at half-filling with disorder as 
obtained by the DMFT with geometric average [48] . Dots represent the numerical solution of the 
DMFT equations. The solid line was obtained analytically by the linearized DMFT. 

4.1 Ground-state phase diagram 

In the half-filled case the ground-state properties of the Falicov-Kimball model are determined 
by states in the center of the band, i.e., at lj = 0. The ground-state phase diagram is shown in 
Fig. [3] In the presence of disorder the metallic phase is seen to be more extended than in the 
pure case. Namely, while for A — the (continuous) Mott-Hubbard MIT occurs at U c — W/2, 
the presence of disorder shifts the transition to larger values of U: the phase boundary between 
the metal and the Mott insulator, Af H {U) in Fig. H lies between < A < 1.70IT and 
W/2 < U < 1.36VF. By contrast, the Anderson MIT, which is also continuous, takes place 
along the phase boundary Af(U) (see Fig. [3]). Again the critical disorder strength increases for 
increasing U and extends between < U < 1.36W and eW/2 < A < 2.03W, where e « 2.718 
is Euler's constant. Obviously the interaction impedes the localization of particles due to the 
scattering at the impurities. For large values of U and A there is no metallic phase. 

The DMFT results shown here confirm that in the disordered Falicov-Kimball model the 
Mott insulator and the Anderson insulator are continuously connected. Hence, by changing U 
and A it is possible to move from one type of insulator to the other without having to cross a 
metallic phase. This scenario is supported by the fact that the Anderson MIT (U — 0) and the 
Mott-Hubbard MIT (A = 0) are not associated with the breaking of a symmetry. 



4.2 Linearized dynamical mean-field theory 

In the vicinity of the MIT the linearized DMFT |52l61j may be employed to calculate the 
transition line analytically. Due to the symmetry properties of A gcom (uj) one finds G(0) = 
— i7rA geom (0) at the band center, i.e., G(0) is purely imaginary. The DMFT self-consistency 
then implies that the hybridization function obeys r?(0) = — inW 2 A geom (0)/16. 
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Iteration of the DMFT equations leads to 



.4 



(n+l) 



(0) 



16 



AW (0)T(e i ) ! 



where 



2 • 



After the geometric average the recursive relations in the linearized DMFT are given by 



(19) 



(20) 



4e^ 3 (0)=4eL(0)^-exp 



While in the metallic phase the quantity A|com(0) increases upon recursion, i.e., Agco-m' (0) > 

^4gcom(0), it decreases in the insulating phase. Therefore, at the boundary between the metallic 

and the insulating solution the recursion does not depend on n, i.e., A^ +1) (0) = A^\0). This 
condition directly determines the MIT transition at A = A(U): 



A/2 



delnT(e) 



-A/2 



(21) 



(n+l), 



1 



w_ 

16 



exp 



i r A ' 2 

- / delnr(e) 

A J -A/2 



w 2 

— exp [I geom (U, A)] 



The integral is evaluated analytically with the result 

2 / a \ 2" 



I seom (U,A) = 2 + In 



-2 In 



(22) 



(23) 



2U 
~A 



arctan 



-In 



A + U 



A-U 



The solution of Eq. (|2"2"j) . shown as a solid line in Fig. [3] is found to agree very well with the 
numerical result. For weak interactions U the critical disorder strength obtained from ((22]) 
increases linearly in U, i.e., A(U) « We/2 + irU/2, since the total bandwidth increases linearly 
with U. For weak disorder A the solution of Eq. ([22]) is given by A{U) « ^JU 2 - (W/2) 2 . This 
agrees with the result obtained from the arithmetic average |48j . which is not surprising since 
for weak disorder both averages must give the same result. 



5 Ferromagnetism and Kondo insulator behavior in the disordered periodic 
Anderson model 

In Sees. Ill and IV we discussed the results obtained for the Mott and Anderson MITs in 
correlated electron systems where the disorder had a continuous distribution, i.e., the probability 
distribution function of the local atomic energies was given by a piecewise continuous function 
of the energy ("box disorder"). However, to model binary-alloy systems a bimodal probability 
distribution function is more appropriate. This type of disorder also leads to a diffusive motion of 
the electrons, and even to Anderson localization, as the disorder with continuously distributed 
energies |50j . But, for sufficiently high disorder strength the binary alloy disorder causes in 
addition a splitting of the band ("alloy band-splitting") in arbitrary dimensions. It is therefore 
important to understand how the opening of two different kinds of gaps - the gap due to the 
alloy band-splitting and the Mott gap - affects the properties of the system. 

In the following we employ the DMFT with arithmetic averaging over the binary alloy 
disorder [55]. As discussed earlier, this approach cannot detect Anderson localization on the 
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level of one-particle quantities such as the density of states. However, for a discussion of the 
effect of band splitting on the correlated electrons such an approach is quite appropriate. In 
any case, the investigation can be extended to include the geometric average [SD] . 

In the case of the one-band Hubbard model with the binary-alloy disorder we found the 
following interesting effects [62163] : (i) for band filling v commensurate with the concentration 
of alloy atoms x, i.e., for v = x or 1 + x, the Mott-Hubbard MIT can take place for non-integer 
values of v, (ii) depending on the ratio between the Mott gap and the alloy gap the correlated 
alloy insulator at non-integer filling can be classified either as a Mott-alloy insulator (when 
the Mott gap is smaller than the alloy gap) or as a charge-transfer alloy insulator (when the 
Mott gap is larger than the alloy gap), (iii) at low filling and for particular values of the alloy 
concentration the Curie temperature increases in comparison with the pure system under the 
same conditions. We now discuss our results for an analogous problem, namely, the effect of 
binary alloy disorder on the properties of correlated / electrons which hybridize with a non- 
interacting band of conduction electrons, as described by the periodic Anderson model (PAM) 

5.1 Periodic Anderson model with binary alloy disorder 

The Hamiltonian of the PAM with binary alloy disorder is given by 

SPAM = %4r C 3<r + Yl (fifth* + e i C L C iJ) 

+ E ( Vc lfi„ + V*flc ta ) + t/ ]T «■ (24) 

ia i 

Here c\ a (c, CT ) and (fia) are creation (annihilation) operators of conduction (c) and localized 
(/) electrons with spin a at a lattice site i. The microscopic parameters in this model are the 
hopping amplitude ty of the c-electrons, the random on-site energies e{ and e|, and V, the 
local hybridization between /- and c-electrons. The Coulomb interaction U acts only between 
/-electrons on the same site. The alloy is modeled by a bimodal probability distribution function 

P(Vi) = xS( yi - y Q ) + (1 - x)S(jh - yo - A v ), (25) 

where j/j = ej?, e{ are independent, random variables with reference values yo = £%, £q- The 
alloy concentration is characterized by the parameter x and the difference between the atomic 
energies of the alloy components by A v = A c , A* , respectively. 



5.2 Alloy-band splitting of non-interacting electrons 

The PAM [Eq. (|24|)] is solved within the DMFT by mapping it onto a corresponding single- 
impurity problem. The arithmetically averaged local Green functions can be written in matrix 
form as 

,ioc^. r„, _ _ ((TrfMftQ)) (T T / CT (r)4(0)) N 
UT rC »/t( )> <T T c»4(0)}, 



Gl oc (T;{y t }) 



They are expressed in terms of local self-energies which appear in the k-integrated Dyson 
equation T, an = Q~^ n — G an . Here Q an is the matrix of local Green functions of the non- 
interacting bath electrons, with 

Q-l _ ( iu n +fi-el V* 
an \ V iu n + [i - e% - Vo 



To understand the effect of the disorder on the physics described by the PAM it is instructive 
to investigate the case (7 = first. For [7 = the corresponding impurity problem is quadratic 
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Fig. 4. Binary alloy splitting of the c- and /-electron density of states for U — 0. a) No disorder 
in the c-electron system (A c = 0). b) Limit of strong disorder, A c — > oo; in this case 2(1 — x)Nl 
non-dispersive /-electron states (dotted line) remain at e' and 2(1 —x)Nl c-electron states are shifted 
to high energies. 



and the functional integrals can be performed analytically. However, in the case of a two-band 
system like the PAM, where /- and c-electrons hybridize, the situation is more complicated than 
in the one-band case discussed earlier, since disorder affects a hybridized two-band system in 
several nontrivial ways. 

We now consider the case where the alloy disorder acts either on the c-electrons or the 
/-electrons, respectively. In the case of c-electron disorder the diagonal local Green functions 
are given by 



Gc 



(G 



cc \-l 
an J 



2rJf 
an 



+ 



1 -X 



(GZ 

V- 7 a 



\v\ 2 gil-A<-- 



G ff = 



1-x 



(GU)- 1 - \V?G? n (GU)-i - (g Jfl 2 _ 4 c 



(26) 



The case of /-electron disorder is obtained by exchanging /f>cin PS)I . Alloy disorder acting 
only on the c-electrons leads to a band splitting of the conduction electrons for large enough 
energy splitting A c . As in the single-band model each alloy subband then contains 2xNl and 
2(1 — x)Nl states, respectively. The c-electron alloy subbands are separated by the energy 
A c . One might expect that, due to the hybridization of c- and /-electrons, a similar effect 
would also occur in the /-electron subsystem. However, this is not the case. Namely, as seen 
from ((26)) a hybridization between the /-states and the 2(1 — x)Nl states from the upper 
alloy c-electron subband is no longer possible for A c — > oo. In this limit a (non-dispersive) 
/-level with 2(1 — x)Nt states appears at the energy e^, in analogy with the case without 
hybridization (V = 0). Consequently, for infinitely strong binary alloy disorder in the c-electron 
system 2(1 — x)Nl /-electron states become localized for arbitrary but finite values of V . So for 
A c — >• oo only 2(1 — x)Nl c-electron states, rather than 4(1 — x) Nl states, are split off from the 
spectrum and are shifted to high energies. We note that, although the band splitting scheme is 
different from the single-band model, the alloy with hybridized c- and /-electrons can still be 
a band insulator for total densities different from integer values (2 or 4). A schematic plot in 
Fig.@]shows the projected density of states, N b (uj) = — ImJ^ £v b (aj)/7r, where b = c or /, for 
a system without [panel (a)] and with [panel (b)] disorder. An analogous analysis of /-electron 
disorder shows that in this case, at large A\ the /-electron band is split into alloy subbands. 
Hybridization between the 2(1 — x)Nl states from the upper alloy /-band and the c-electrons 
is again prevented when A? — > oo. Therefore, the corresponding fraction of the c-electron band 
is unchanged, i.e., remains at the same energies as in the non-disordered case. We thus see that 
even in the absence of interactions binary alloy disorder affects a hybridized two-band system 
and a single-band system in quite different ways. 



Will be inserted by the editor 



13 



0.3 
0.2 
0.1 


0.2 
0.1 


0.2 
0.1 


0.2 
0.1 







J /X A c = 0.0 




'/^"^\^^^^ A c = 2.0 . 


A c = 4.0 




-2 


2 4 6 8 10 
CO- (J, 



Fig. 5. Spectral function of c-electrons in the PAM for different A c at x — 0.3 (other parameters: 
U = 1.5, V = 0.5, £q — e o = 3-25 and n to t = 1.3) obtained within QMC and maximal entropy at 
T — 1/60. By increasing A c a pseudogap opens, which becomes a real gap for T — > 0; after Ref. |64| . 



5.3 Kondo-insulator of interacting electrons at fractional filling 



In analogy to our findings for the Hubbard model with binary alloy disorder 62 63 an opening 
of a correlation gap, obtained by increasing the alloy band splitting A c , is found for the PAM 
[64] ; see Fig. [5] This is caused by the splitting of the c-electron band due to binary alloy disorder 
and the correlations between the /-electrons. Namely, when the energy splitting A c is much 
larger than the width of the c-electron band the total number of available low-energy states is 
reduced from ANl to [4 — 2(1 — x)]Nl — 2(1 + x)Nl, whereby the filling effectively increases 
by a factor of 4/ [2(1 + a;)], such that n^ t = 2n tot /(l + x), if n tot < 2(1 + x). For the filling 
"-tot = 1-3 studied in Fig. [SJ the concentration x — 0.3 is a special case since then n^ t = 2. 
The system is then effectively at half- filling and behaves as a Kondo insulator at large U, 
A c , and low temperatures. The transition from a metal to a Kondo insulator at non-integer 
filling predicted here for the PAM [64] is a counterpart to the Mott-Hubbard metal-insulator 
transition at non-integral filling in the one-band Hubbard model discussed in [62163] . 



5.4 Disorder induced enhancement of the Curie temperature 



It is well known that itinerant ferromagnetism occurs in the non-disordered Hubbard model 
(fT]) only off half- filling provided the DOS is asymmetric and peaked at the lower edge [65166] , 
While the Curie temperature increases with the strength of the electron interaction one would 
expect it to be lowered by disorder. However, our investigations show that in some cases the 
Curie temperature can actually be increased by binary alloy disorder [62167] . 

In the case of the PAM with binary alloy disorder we also found that under certain conditions 
the Curie temperature can actually be enhanced [63]- Indeed, as shown in the right panel of 
Fig. [6l the Curie temperature for the transition to the ferromagnetic state in the PAM is a non- 
monotonic function of the alloy concentration x. In particular, the behavior is quite different 
for disorder acting on the /- or the c-electrons. 
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Fig. 6. Curie temperature in the PAM as a function of alloy concentration x and energy splitting 
(left column) and A c (right column) for n to t = 1.3 and e§ ~ e o = 3.25. Strong c-electron disorder 
enhances T c compared to its values at x — or 1; after Ref. [64] . 



5.4.1 /-electron disorder 



As found by Meyer |68| the presence of /-electron disorder always reduces the Curie temperature 
relative to its non-disordered values at x = or 1. For strong enough disorder T c eventually 
vanishes, e.g., at x = 0.28 and x = 0.75, respectively, for A' = 1.7 (right panel, left column 
of Fig. IS]) . This is due to the splitting of the /-electron band at large A? which increases the 
double occupation of the lower alloy subband; this reduces the local moment of the /-electrons 
and thereby T c . 



5.4.2 c-electron disorder 



Interestingly, c-electron disorder leads to a much more subtle dependence of T c on concentration 
x. Namely, for increasing energy splitting A c there are, in general, three different features 
observed, the physical origin of which will be discussed in more detail later: 

(i) at x — 1, i.e., in the non-disordered case, T c is reduced, 

(ii) a minimum develops in T c at x — n to t — 1 > 0; 

(iii) T c is enhanced over its non-disordered values at x — or 1. Altogether this leads to a global 
maximum in T c vs. x. While the decrease of T c at x = 1 is a simple consequence of the reduction 
of the energy difference between the /-level and the c-electron band, e c — — e§ — £q — A c , 
for increasing A c , the latter effects are more subtle. 
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5.4.3 Origin of the maximum in T c 

We will now explain the maximum in T c vs. x. It can be understood within the following model 
based on an ansatz for the Curie temperature, T C (U, V, fi) = T®{U, V, ^)F c (fi — eg)F^(^ — £q), 
which implies that the formation of local /-electron moments (F*) is assumed to be independent 
from the c-electron mediated ordering of those moments (F c ). In fact, for the RKKY model 
this ansatz can be microscopically justified within a static mean-field theory. The two functions 
F c , Ff are determined by T c calculated within DMFT for the non-disorder case at fixed fj, — £q 

or /it — £g, respectively; they are shown in Fig. [Tfa) and[7]Jb) for one set of parameters. The 
prefactor Tj? is determined by the requirement that the dimensionless functions F? and F c be 
equal to one at their maxima. We note that F^ (// — Eq ) has a maximum when the /-level is 
half-filled (/z = £q + U/2), i.e., when the local moment is maximal. 

The Curie temperature in the presence of c-electron disorder can be estimated by averaging 
over the c-electron part, F c , giving rise to the disorder-dependent function F c (x, fi — Sq) = 
[xF c (fi — £§ + A c ) + (1 — x)F c (fi — £q)]. The linear dependence on the alloy concentration 
can again be justified microscopically within a static mean-field theory for the RKKY model, 
where T c depends linearly on the DOS at the chemical potential. T c is now determined for 
each concentration x. We calculate [i, which is an implicit function of x, in the non-hybridized 
limit (V = 0) within a rigid band approximation. The dependence of the resulting functions 
T c {x 1 fi — £q) and F^(/j, — Eg) on x are shown in Fig. (7[c) for A c = 2.0. In general F? (// — Eg) 
has a global maximum at those values of x for which the /-level is half-filled [see Fig. E{c)]. 
By contrast, F c (x, /i — eg) is characterized by a wide minimum, related to the formation of the 
pseudo-gap in the interacting DOS seen in Fig. [S] This minimum reaches zero, i.e., J- C (x,fi — 
£§) = 0, for a finite range of x values as shown in Fig. \T[c). The resulting T c (x) obtained by 
the product of these two functions agrees remarkably well with the numerical result obtained 
by DMFT as shown in Fig.EJd). 

5.5 Summary 

In this paper we showed that the interplay between disorder and many-body correlations 
can lead to quite unexpected, often counter-intuitive, behavior. We first explored the zero- 
temperature phases of the Anderson-Hubbard model within the dynamical mean-field theory 
in combination with a geometrical average over the disorder. This allows for a unified descrip- 
tion of Anderson and Mott-localization on the basis of one-particle correlation functions. The 
paramagnetic phase diagram shows reentrant metal-insulator transitions caused by the interac- 
tion and disorder, and the Anderson and Mott insulating phases were found to be continuously 
connected. In the presence of antiferromagnetism a stabilizing effect of the simultaneous pres- 
ence of interaction and disorder was discovered, which leads to a new antiferromagnetic metallic 
phase. An overall similar behavior concerning Anderson and Mott insulating phases is found in 
the Falicov-Kimball model. It is expected that the very interesting parameter regime of strong 
interactions and strong disorder, which is not easily accessible in correlated electron materials, 
can also be realized with cold atoms in optical lattices. 

In a second application we investigated ferromagnetism and Kondo insulator behavior in 
the periodic Anderson model (PAM) in the presence of binary alloy disorder. Far away from 
half-filling the PAM shows ferromagnetic order. For strong enough binary alloy disorder the con- 
duction band splits and the correlations among the /-electrons lead to a non-trivial dependence 
of the Curie temperature on the alloy concentration. Upon decreasing the alloy concentration 
the local moments increase which raises T c , but at the same time the opening of a gap in the 
alloy Kondo insulator at non-integral filling leads to a decrease of T c . In effect this causes the 
Curie temperature to behave non-monotonically as a function of the alloy concentration, with 
a global maximum in T c which can be drastically larger than in the absence of disorder. The 
effect is predicted to occur in /-electron materials with alloy disorder in the conduction band, 
and also in ultracold fermionic atoms in optical lattices trapped by harmonic potentials in the 
presence of random binary disorder. 
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X X 

Fig. 7. (a) F c (n — eg), (b) F' (fx — Eg) appearing in the ansatz for T c in the disordered PAM (see 
text) calculated for A c — 0. (c) J-" c (/x — s§) and F^(ii—£q) for Z\ c = 2.0; other parameters as in Fig. [6] 
(d) Comparison of T c obtained from the ansatz and within DMFT; after Ref. |64| . 
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